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ABSTRACT 


The effect of mobility ratio on areal sweep efficiency and con- 
ductance ratio has been made for direct and staggered line drive well 
patterns. This study was made using miscible fluids, artificially 
consolidated sandstone models, and a photographic method of data 
recording. The mobility ratio was varied between 0.1 and 12, which 
is the range common to water flood operations. 

A method of predicting the performance of waterfloods in 
stratified reservoirs using conductance ratios is treated. This 
technique is not limited in number of strata and indicates one method 


of considering a free gas saturation. 
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INTRODUCTION 


The importance of secondary recovery techniques in petroleum 
engineering is continually increasing and is taking an increasing per- 
centage of reservoir research activity. Since primary reservoir deple- 
tion normally will leave approximately 75 per cent of the original oil 
in place, the reason for this increasing role of secondary recovery is 
obvious. 

Perhaps the oldest and most used method of recovering oil by 
secondary techniques is the waterflood. This involves injection of 
water into the reservoir to displace the oil. Water injection requires 
expenditure of energy and funds, therefore, reliable methods for pre- 
dicting waterflood results are necessary. 

Extensive study has gone into conformance factor or sweep 
efficiency predictions for various well patterns at varying mobility 
ratios. This permits predictions for single stratum reservoirs. 

Most reservoirs are composed of several strata, however, and 
performance predictions for stratified reservoirs are not feasible 
with only sweep efficiency data available. In predictions for stratified 
systems, it is necessary to place performance ona time basis. This 
requires that flow potential or flow rate for a pressure drop be avail- 


able for the system under investigation. 


Flow potential data (conductance ratios) are available for all 
(AAUR Se ; ; 
regular patterns with the exception of the line drive. The pur- 
pose of this study was to obtain line drive conductance ratios and to test 
methods for approximating conductance ratios analytically. A method 


of predicting performance of a stratified system using these data is 


also presented. 


WATER FLOOD PREDICTIONS 


In studies of petroleum reservoirs, much qualitative and quanti- 
tative use has been made of models to more completely understand the 
phenomenon of oil production. Since no two reservoirs are exactly the 
same, model studies are limited in scope. Therefore, ina rigorous 
interpretation, models can only be used for studies of a particular 
field; however, by generalizing parameters, study of a reservoir seg- 
ment will yield data for qualitative application in a number of situations. 

A great deal of work has gone into the determination of scaling 
criteria for model studies. ole Co ame Probably the most 
complete analysis covering general situations has been presented by 
Geertsma, Croes, and Schwarz. ot A strict interpretation of this 
analysis, however, leads one to the concl usion that the model must 
be identical to the prototype in every respect. Since this is completely 
impractical, it is necessary to decide which reservoir factors are of 
greatest importance and construct the model accordingly. A detailed 
description and application of model scaling for edgewater drive reser- 
voirs has been presented by Caudle. 

The efficiency of a secondary recovery project, in which a fluid 
(usually water) is injected to displace the oil in the reservoir, can be 
defined by the volumetric sweep efficiency. In reservoir studies, 


models are normally two dimensional representations and the displace- 


ment efficiency may be defined in terms of areal sweep efficiency. This 
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is the ratio of the area swept by the injected fluid to the total model area. 

Generally, the areal sweep efficiency is determined by the well 
pattern geometry and the mobility ratio. This assumes that the reser- 
voir is homogeneous and that gravitational forces are not present. The 
mobility ratio is defined as the ratio of the sum of the mobilities of the 
fluids behind the front to the sum of the mobilities of the fluids ahead 
of the front. The mobility (}\) of a fluid in a porous media is expressed 
as the effective rock permeability (k) to the fluid of interest at its 


saturation divided by the fluid viscosity ({). 


iN 
(k/u) swept S ae 
M = 3 = = = Mobility Ratio (1) 
(k/u) unswept LY 


Early model studies utilized the electrolytic model to represent 
waterfloods. This model uses the principle of ionic flow to simulate 
fluid flow. Use of this model is limited to mobility ratios of one. 
Muskat, Wyckoff, and een a were among the first to utilize this 
method to determine breakthrough sweep efficiencies for various well 
patterns. Ina later publication, Muskat and Wa contee analyzed the 
previous electrolytic model poauiese | and concluded that the shape of 
the flood front and the breakthrough areal sweep efficiency were depend- 
ent on the geometry of the system and not on pressure differential, well 


2a, Co. coer 
Spacing, or rock properties. It is noted that Muskat indicates 





that the results obtained in the electrolytic models do not consider the 
effects of two fluids with different viscosities. In other words, the elec- 
trolytic model is valid for a mobility ratio of one only. 

The electrolytic model and the similar but more sophisticated 
potentiometric model have since been used in the study of various reser- 
voir conditions. Limited use has also been made of heat flow models to 


oe 
simulate fluid flow with a mobility ratio of one. 


Sweep Efficiency 

Early studies on the effect of mobility ratio on areal sweep effi- 
ciency utilized potentiometric models. Avoneisky” reported results of 
a potentiometric model study of mobility ratio effects. A stepwise pro- 
cedure was utilized to vary the ratio of the depth of the electrolyte 
behind the front to that ahead of the front. This variation in depth made 
the electrical conductivity ahead of the front different from that behind 
the front, thereby simulating fluids of different mobilities flowing ina 
porous media. After computing a new frontal position, the electrolyte 
depth was changed to fit the new frontal position and the process was 
repeated. In thisiterative manner, sweep efficiencies for mobility 
ratios other than one wereapproximated by advancing the front ina 
series of steady-state steps. 

The effect of mobility ratio on areal sweep efficiency has also 


Z 
been studied with an electric analog model by Nobles and Janzen. 
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6 
This study utilized a resistance network and an iterative process similar 
to eewotsky's. to simulate frontal movement. 

Cheek and Weazie. studied the effects of mobility ratio on areal 
sweep through the use of a fluid mapper. This model uses fluid flow 
between two closely spaced parallel plates to simulate fluid flow ina 
porous media. Results are reported for five spot and direct line drive 
well patterns. 

Slobod and Cerealle introduced the X-ray shadowgraph technique 
to study fluid flow phenomenon. This new technique used a porous plate 
and fluid flow with one of the fluids containing an X-ray absorbing mate- 
rial so that a contrast between displaced and displacing fluid could be 
observed on X-ray shadowgraphs taken during the experiment. By 
taking successive X-ray exposures, a complete history of the flood was 
obtained. 

Dyes, Caudle, and peeves used the X-ray shadowgraph 
method to illustrate that a significant amount of production could take 
place after initial breakthrough. The results of this study demonstrated 
that, with most water floods, an ultimate areal sweep efficiency of 85- 
100 per cent could be obtained for various mobility ratios at an economic 
water-oil ratio. The only previous reports which indicated consideration 
of production after breakthrough were made by Eee oe for a mobility 
ratio of one and Fay and eats ae for a mobility ratio of 4.0 using 


numerical techniques. 







ete eG ei, 
ees Sey 
Sa) os hi < qs ~ pe 


222 2S) GQ =e 2 © eo —s 
© ——T ee = =e a ae 


eh oo = ee an 
~~ Oe ee = = & —- 
—_ 2 a h 6«8@ml ee => ~,,. =, 
= —_—_> & ea a= ——= 
Ss © — 


The X-ray shadowgraph has been used extensively to study a 
variety of reservoir conditions at varying mobility ratios. 7 Os Tt ee 
Included among these studies are five-spot pattern floods, stratified 
five-spot floods, and the effect of fractures on flood performance. 
reeeermann © utilized a method similar to the X-ray shadow- 
graph for studying fluid flow. This method consisted of artificially 
consolidated sand models covered with clear plastic. One of the fluids, 
either displaced or displacing, contained a dye which permitted visual 


observation of the front as it progressed. Photographs taken during 


the flood permitted calculation of areal sweep efficiency. 


Conductance Ratio 

The use of areal sweep efficiency data permits ultimate re- 
covery predictions for single stratum reservoirs. In order to predict 
performance of stratified reservoirs and also to predict recovery of 
single stratum reservoirs ona time basis, it is necessary to know the 
change in the flow potential or flow rate per unit pressure drop for the 
system under study. 

Prediction of flow potential was first discussed by mueieres in 
terms of conductivity (production rate per well per unit-pressure 
differential). This permitted predictions for systems with mobility 
ratios of one. 


2, 
Aronofsky and Ramey first introduced the concept of conductance 





ratio ina potentiometric study of mobility ratio effects ona five-spot 
water flood. The conductance ratio (y) is defined as the flow rate (q) 
per unit pressure drop (AP) between the injection and production wells 
at a specified time to that conductivity with only the original fluid pre- 


sent. 


os q/AP) time t = Conductance Ratio (2) 


(q/AP) initial 


Analytical methods of studying the change in flow potential has 

been presented in terms of dimensionless injectivity by various methods. 
30 ar: _ 
Prats, Matthews, Jewett, and Baker have computed injectivity for 
five-spot systems by considering that the reservoir consists of several 
ee I 
layers, each having a uniform but different permeability. Hauber 
presented another approximate method which considers that flow uti- 
lizes the stream tubes and isc-potential lines obtained by the solution 
of differential equations for a mobility ratio of one. 
2 mee 

Deppe presented a method to approximate conductivity for vary- 

ing mobility ratios and for varying well patterns. In the five-spot 
iZ ; — 

pattern used as an example, Deppe essentially divided the pattern 
into two regions of radial fiow. The first region was around the injec- 
tion well and covered almost one-half of the area. The second radial 


flow region was the same except that it was around the production well. 








7 

The results of Denne study provided injectivity (conductance 
ratio) curves which were not too unreasonable for five-spot and nine- 
spot well patterns. Closer agreement was obtained for unfavorable 
mobility ratios, which is normally true when using this type of approxi- 
mation methods. 

In an attempt to improve on the work of Deppe, this study has 
divided the pattern into 3 regions as illustrated by Figure 1 fora 
staggered line drive pattern. The first region surrounds the injection 
well where the largest pressure drop is encountered. This region is 
considered to have radial flow and has a radius equal to one-fourth the 
well row spacing. The second region is considered to be linear, with 
only slight pressure change, and extends to the areal sweep at break- 
through. The third region is one of radial flow around the production 
well where the pressure drop is considered to be a function of frac- 
tional flow from the swept region. 

Conductance ratio for the system is computed by considering 


a constant flow rate which reduces the equation to 


. = 
eee (3) 


AP. 





Pressure drop across the system is obtained by summing the pressure 


drops in the three regions at the time of interest. 
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Production Well 


Linear Flow 
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Injection Well 


et Radius of radial flow area surrounding injection well 
rp Radius of radial flow area surrounding producing well 
fo! Angle open to flow from swept region after breakthrough 


FIGURE 1. STAGGERED LINE DRIVE -- ASSUMED FLOW REGIONS 





)] 

With only single fluid flow in the system (or at initial conditions) 

the pressure drop in each region can be represented by modifications of 
Darcy's Law as 


Gola) ay 4 
A” aeons (4) 


(5) 


id 
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_ qin rp/rw 
Ser a Se AU (6) 
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which when added, gives an initial pressure drop for the system of 


AP. xa on ling ri/rw + ye 28 L/w + In rp/rw) C7) 


After the flood has started and two fluids are flowing in the model, the 
pressure drop can be represented as a summation of the pressure drop 
in each fluid region. When the front is in the first region the pressure 


drop therein can be represented as 


AP _ qin rf/rw + 4 Tray Tt 
At — eeOeh XT. 08h h 


(8) 


_ q (1/M) 1n rf/rw + In ri/rf 
- ~ AT. 08h 
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When the front is in the second region the pressure drop therein is 


_ qs q(L- 4) 
APpt = X,1.127aw * \ 1.12 7hw 


(9) 
<r 1 s/w t+ L/w 
q ela 


After breakthrough, the flow rate in the third region must be summed 
for parallel flow since two fluids are flowing into the producing well. 
Consider that portion of the well illustrated in Figure 1 or one-quarter 
of the well. Let a equal the angle open to flow from the swept region 
and (1/2-a) equal the angle open to flow from the unswept region. 


Then the fractional flow from the swept region is 


rot 


= Ss 
fs = ycaph (i/ 2a) ie 
Ss u 
then 


ae ISTRET uh) 





13 


The pressure drop in the swept and unswept region is 


fs ¢@ Inge rue / Zz 


APCs ~ ; X_ a (12) 
Ss 





(l-fs) q ln rp/rw 7/2 


em TT UW. Gal Soe bX (n/2-a 


(13) 


| For parallel flow, the pressure drop given by equation (12) must equal 





that of equation (13). Therefore, the pressure drop in the third region 
may be obtained from either the equation for the swept or the unswept 


portion. By substituting equation (11) into (13) we obtain 


(l-fs) q ln rp/rw 


“Pot ~ 7.08h X, [l-fs/[M-fs(M-1)] | 


(14) 


In summing equations (8), (9), and (14) the following conditions 


apply: 
1, When Es<Es ,: O/w = (0, fs = 05 and rt = i( hs). 
2. When Es, <Es<Es,,: rt = ri, 1s = 10> and 4/ Ww = s). 
3. When EsvEs_._¥§: rf= ri, 2/w= L/w, and fs = f(Es). 
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The total pressure drop at any time is then 





AP, = 17.08 asm In rf/rw + ln ri/rf 


+ 


6. 28 Ess : 1] he CeO ee (15) 


=F 


(l-fs) In rp/rw 
T-571M-is(M-1)] 


The conductance ratio from equation (3) becomes the ratio of equation 


(7) to equation (15) or 


yaa AP./AP, = (In ri/rw + 6.28 L/w + 1n rp/rw)/ 


cum) In rf/rw + 1n ri/rf + 6. 28 {xr - y L/w (16) 


(l-fs) ln rp/rw 


+ 6,28 L/w + sie eyaIiets| oan 


Values for L/w were obtained by using both arbitrary values and by 
assuming that the pressure drop across a linear system is equal to 

that across a line drive with well radius equalto ri. Muskat's equation 
for flow rate in a line drive is: 


a d/a - 1.17 + 0.637 log a/ri (17) 
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Setting equation (17) equal to (5) and solving for L/w yields 


Co Apa ntti oni togalr (18) 


which is valid when d/a>l. Calculations were performed by programming 
the above for the Controlled Data Corporation 1604 digital computer. 
Appendix D contains the flow diagram and FORTRAN language program 
utilized in these computations. A conversion table for the FORTRAN 
symbols is also provided in Appendix D. 

In the calculations for conductance ratios, various values of L/w 
were assumed in an attempt to match the model study results. Using 
equation (18), L/w was 0.357. Other values used for L/w were 0 and 
0.1. As can be seen from the plots of Appendix D, varying L/w made 
only small changes in the conductance ratio. Comparison plots for L/w 
of 0.1 and as obtained from the model studies, indicate results similar 
to those obtained by Deppe. Agreement is fair for unfavorable mobility 
ratios and not too good for favorable mobility ratios. 

This method can be used for rough approximation of conductance 
ratios for varying well patterns by modifying the initial values in the 
program to correspond to those of the pattern under study. However, 
it is believed that fluid flow model study results would be more reliable. 

The effect of mobility ratio on conductance ratio has been 


studied for various well patterns using fluid flow models. Caudle and 
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16 
weed have published data for the five-spot, the inverted skewed seven- 
spot has been studied by Hickman, a8 and the inverted nine-spot was 
studied by Watson. Zi This study provides data for line drive well 


patterns. 


Stratified Reservoirs 

In a reservoir composed of strata of varying permeabilities 
and porosities, performance prediction becomes more complex than 
for a single stratum. To illustrate, consider a two strata system of 
equal characteristics with fluids of equal viscosity and with the perm- 
eability in the first stratum twice that of the second. From Darcy's 


Law 


AP 
a (19) 


< 
Tt 
| 
cin 


it can readily be seen that the first stratum will have a frontal advance 
twice that of the second strata and that at breakthrough for the first 
stratum the second stratum will have a volumetric sweep of only half 
of the breakthrough sweep efficiency. In this simple example, a mis- 
cible system was assumed in order that mobility ratios of one could 

be used for each strata. For mobility ratios other than one, the linear 
relationship of flow velocity to permeability no longer exists ina 
stratified system and then becomes a function of distance traveled or 


area swept. 
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In order to apply the data obtained in this study to a stratified 
reservoir, it is necessary to put the performance of each strata on an 
equal time basis. This is most easily done through the ry the con- 
ductance ratios obtained in this and other studies. It is noted that this 
method is developed and accurate for separated strata or for strata 
without cross flow. When cross flow between strata exists and the 
mobility ratio is greater than one, this method provides a good 
approximation for predicting water flood results. A similar method 
using theoretical dimensionless injectivities has been presented by 
Pratts, etal. 0 The method developed herein permits performance 
prediction as accurately but more simply than the method described 
by Prats, etal. ze 
To develop this time relationshlp, consider Darcy's Law (19) 
and let 
v = q/A 
and 
A = Vd/¢£ 


Darcy's Law can then be written in terms of conductance as 


q/AP = -——— = conductance (20) 


Consider now two adjacent parallel strata in a reservoir with single 


fluid flowing; the variables between these two strata are permeability (k) 





and displaceable volume (Vd). 
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Therefore the conductance between a stratum 1 and any stratum 


j will vary as their permeability and displaceable volume ratios vary 


or 


ha ae 
BP; Ri Vd) OP, 


By definition, 


y = (q/AP), / (a/4P), 





then 
Ps ed - q 
|a,/aP, = ¥; vd = 
J Ue J Kk) V4, LAP) 4; 
and 
Seiya i ai | 
emcee ae Sich abysll 


solving for dt, integrating, and letting c equal AP, | 


k, Vd, 


Citas KVa, nu dwi 


(21) 


(22) 


(23) 


(24) 


(25) 








Lo 
Other equations needed to express stratified reservoir perform- 
ance are for total flow from the swept region (fs) and total conductance 
hat1o (+). 


By definition, 


fies = yee 26 
3 eC ae oy Ge (26) 
then 
s. = fs. q. Zi 
a ; 4 (27) 
and 
AP, 
=, sale 28 
EE GE Le [aPT, (28) 
from (7) 
eq: 
ia. | = Kk, va, la, 1; (29) 
then by substituting into (26) we have 
yah Ne DE Jee ol 
Ee _ 
Ae J 


The expression for overall conductance ratio can be derived in this 


Same manner as 
Seley ch 
Ne eT eat 


Y > Se va. 
ingee 


(31) 


An example performance prediction for a stratified reservoir 
using this method is described in Appendix E. A possible method of 
considering an initial free gas saturation is also mentioned in Appendix 


E. 





EXPERIMENTAL APPARATUS 


Model Theory 

The majority of oil reservoirs are drilled with a regular pattern 
of wells. Probably the most common pattern utilized in this country is 
a square array of wells ona 40 acre spacing. This well pattern, illus- 
trated by Figure 2, readily lends itself to secondary recovery operations 
normally referred to as the Direct Line Drive, which is obtained by con- 
verting alternate rows of wells to injection wells. 

The Staggered Line Drive is also obtained by converting alternate 
rows of wells to injection wells; however, the wells are not drilled ina 
square array but rather inan offset fashion. The well pattern for the 
Staggered Line Drive system is shown by Figure 3. 

With the regular well systems shown by Figures 2 and 3, it can 
readily be seen that there will exist numerous mirrow images divided 
by lines of symmetry. These lines or planes of symmetry represent 
invarient streamlines or boundaries across which no flow will occur, 
as long as production and injection quantities are uniform. J Since the 
planes of symmetry are in effect flow boundaries in a homogeneous 
system, it is apparent that only the smallest segment bounded by 
symmetry planes need be considered. Figures 4 and 5 show enlarged 
views of the model and its relationship to the line drive systems studied. 
As can be seen, the model utilized represents a line drive pattern where- 


in the row spacing equals the well spacing. 
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Model Construction 

The model utilized in this study was an artificially consolidated 
sand plate constructed generally as described by Weteon and Hickman. ; 

The artificial sandstone was composed of a clear epoxy resin and 
60-100 mesh Ottawa sand. Actual construction of the artificial sandstone 
has been reported in detail by Watson, "4 Hickman, ee and Caudle. ? In 
order to avoid any possible loose streaks which may have occurred during 
packing of the artificial sand a minimum of 2 inches was cut from the 
model sides by using a diamond saw. To seal the model, an extremely 
viscous epoxy resin was spread over all exposed sand faces. The high 
viscosity of the epoxy resin was necessary to prevent imbibition of the 
coating into the model and was obtained by applying the coating just 
prior to it being completely "'set-up."' 

Injection and producing "wells'' were constructed from 1/2 inch 
thick acrylic plastic plate as shown by Figure 6. Swagelok fittings for 
1/8 inch tubing were utilized to connect the "wells" to the injection 
systems. Prior to attaching the plastic plates to the model, well bores 
were reamed out at the well corners. 

As seen from Figure 6, the model utilized was 32 cm in length, 
16 cm in width, 0.6 cm thick, and had a well bore diameter of 0.2 cm. 


This represents a d/rw ratio of 320. 
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The fluids utilized in this study were miscible in all proportions. 
The use of miscible fluids to simulate reservoir fluids in model studies 
of mobility ratio effects has been well established with experimental data 
by Slobod and eamae and in other publications. Bal 14133 

By using miscible fluids in this study, there are no effects of 
capillary forces as will be found in immiscible systems. The absence 
of capillary effects eliminates the problem of saturation gradients in 
determining the effective mobility ratio. Also by utilizing miscible 
fluids, the effective permeability to each fluid is the same, thereby 
eliminating the need to determine relative permeabilities. This system 
therefore models a water flood with little or no flow of oil behind the 
Pree ne. Srieiwobility ratio, GSH) SYeP nen reduees to the 


(k/u) unswept ’ 
viscosity ratio of AE RURISIE Pe . Therefore, only the viscosity of the 


(Uu) swept 
fluids needed to be changed in order to vary the mobility ratio of the 
system. 

To obtain fluids of the desired viscosities, a mixture of gly- 
cerine and water was used. The quantity of glycerine added for the 
desired viscosity was determined on a weight per cent basis through 
the use of standard chemical tables and checked with an Ostwald 


viscosimeter. 


In the photographic technique used in this study, one of the fluids 
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contained a dye to establish a contrast between the displacing and the 
displaced fluids. Since the epoxy used to consolidate the sand renders 
the grains oil wet, oil soluble dyes could not be used. It was found that 
an ordinary household food coloring used in about 5 per cent volumetric 


proportions with the water obtained the desired contrast. 


Injection System 


The constant rate pumps, utilized and described by ricer 
and mateen” were used in this study. This pump is basically a hydraulic 
cylinder using a screw type mechanism to move the piston. It is driven 
by a gear and chain arrangement from an electric motor through a 
variable speed reduction gear. By changing the gears in the drive mech- 
anism and/or the speed output of the reduction gear, the desired injection 
rate could be obtained. 

A pressure gauge was located at the injection well. The producing 
well was permitted to discharge at atmospheric pressure, thereby per- 


mitting the gauge to indicate the pressure drop across the model. 


Data Re cording 


A 35 mm camera was used to record the frontal position, pressure 
reading, and time during the runs. The camera was actuated by a D.C. 
timing motor which advanced the film and made one exposure per minute. 
Exposures could be made at longer intervals by using a ''sign flasher'"' to 


interrupt current flow to the timing motor. This was done during the 
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slower runs. A red filter was used to obtain better contrast between 
the blue and clear fluids on black and white medium speed panchro- 
matic film. A more complete description of the photographic and 
developing technique has been previously written by Caudle. : 

Figure 7 is a schematic diagram of the equipment used in this 
study. 

Typical photographs are shown in Figures 8 and 9. Figure 8 
is an exposure of a direct line drive at a mobility ratio of one and Figure 
9 is for a staggered line drive at a mobility ratio of twelve. Both expo- 


sures indicate frontal positions shortly before breakthrough. 
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PHOTOGRAPH OF MODEL DURING A DIRECT LINE 
DRIVE RUN; .M = 1. 
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PIGURE 9. 


PHOTOGRAPH OF MODEL DURING A STAGGERED 
LINE DRIVE RUN, M = ld. 
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EXPERIMENTAL PROCEDURE 


The performance of both direct and staggered line drive patterns 
were studied for mobility ratios between0. land 12 to determine the 
effect on areal sweep efficiency, fractional flow, and conductance ratio. 

The model was initially evacuated by means of a vacuum pump. 
After complete evacuation, the model was saturated with water. At 
this time, all lines and the pressure gauge were connected to the model 
and the injection pump (Figure 7). The injection pump was loaded with 
dyed water as the initial displacing fluid, the model, pressure gauge, 
and clock were placed ona light table beneath the 35 mm camera and 
the direct line drive runs were started. 

As injected fluid moved into the model, the automatic camera 
recorded the advancing front positions and injection pressures at regular 
intervals. The camera was stopped after approximately 95 per cent of 
the model had been swept. In either case, injection was continued until 
complete displacement of the ''oil analog" had occurred. 

In the second run, the injection system was loaded with fluid of 
higher viscosity and the above process was repeated with the clear 
viscous fluid displacing the dyed fluid. By alternating injection between 
clear fluid of varying viscosity and the dyed fluid, the range of mobility 
ratios from 0.1 to 12.0 was covered in five runs. 


After completion of a set of runs covering the range of mobility 


oN) 
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ratios under study, the roles of the injection and producing wells were 
reversed and the five runs were repeated. This procedure was then 
followed for the staggered line drive system. The two runs at the 
same conditions were made in opposite directions to average out any 
difference in the pressure drop around the individual wells. 

The flow rate was set prior to the start of each run and was 
maintained constant throughout the run. The rate was chosen so that 
the injection pressure would not exceed the gauge scale and so thata 
sufficient number of photographic exposures would be made during the 
run to permit interpretation. The areal sweep efficiency is independent 
of flow rate in the line drive floods. Since the conductance ratio 
(Equation 2) is a dimensionless ratio, the actual rate utilized is of no 
concern. It is important, however, to ensure that the flow rate remains 
constant during each run so that the time scale can later be converted to 
pore volumes of water injected and the conductance ratio computed on 


the ratio of pressure drop. 





ANALYSIS OF DATA 


Areal Sweep Efficiency 

All data was recorded on film during each run. After developing 
the film, a microfilm reader was used to enlarge the recorded image. 
Tracings were then made of the model and the frontal positions during 
sweepout. Appendix A contains typical tracings of frontal positions at 
various mobility ratios for both the direct line drive and the staggered 
line drive. The area of the model and the area behind the frontal posi- 
tion were obtained with a polar planimeter. The areal sweep efficiency 
was then obtained for each frontal position by dividing the area swept by 
the area of the model. Tables of data are contained in Appendix B. 

Areal sweep efficiencies were plotted versus an elapsed time 
scale. Since runs at the same mobility ratio, but in opposite directions, 
were not necessarily made at the same rate, a time correction factor 
was obtained by plotting areal sweep data for the second run on the 
first curve obtained and then noting the time correction necessary which 
would make the curves most closely coincide. This permitted easy 
correlation of data for runs of the same mobility ratio made in opposite 
directions. 

Since each run was made at a constant rate, it is easily seen that 
the time scale of the basic data plots is directly related to the volume of 
displacing fluid injected by a constant ratio for each particular rate. 
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Prior to breakthrough of injected fluid, a plot of areal sweep versus 
volume injected is a straight line. By noting the area swept and the 
elapsed time at breakthrough, all subsequent time values may be con- 
verted to displaceable volumes of injected fluid by the ratio of area 
swept to elapsed time at breakthrough. Appendix B contains plots of 
all mobility ratios studied for the two well patterns considered, with 
tables of data derived from these plots. 

Plots of areal sweep versus displaceable volumes injected can 
be found in Appendix C. The displaceable volume injected is a more 
useful term than actual volume (or time at constant injection rate) due 
to its dimensionless nature. Dyes, etal” has defined displaceable 
volume as the product of the hydrocarbon pore volume of the unit 
pattern and the volumetric displacement efficiency. 

It is noted that areal sweep data from the first 5 runs did not 
correspond to other runs for the direct line drive made in the opposite 
direction as regards areal sweep efficiency. After study of the photo- 
graphed flood front shapes, it was noted that an air pocket had been 
present near the production well during these runs, thereby distorting 
streamlines and causing invalid areal sweep data. This error was 
noted after completion of all runs and an attempt was then made to re- 
make these runs. However, the injection well became "plugged" during 


these ''reruns'' and areal sweep data for the direct line drive for mobility 





——————— 


on 
ratios of 3.13, 0.106, and 12.1 is based on one run only. The cause 
of this 'well plugging" is attributed to impurities in the injected fluid, 
which could have been prevented by installing a filter on the injection 
line. This air pocket did not interfere with runs subsequent to run 5 
since a number of days elapsed between run 5 and 6 and the model was 


routinely evacuated before making run 6. 


Conductance Ratio 

The conductance ratio is defined by Aronofsky and Romney 
@/aP), / (a/P), which reduces to AP, | AP, when the flow rate is 
maintained constant. Pressure data taken from the injection well pro- 
vides AP ata particular areal sweep efficiency and permits direct 
calculation of the conductance ratio. 

From the above equation, it can be seen that if the pressure 
increases during a run, the conductance ratio will decrease and a de- 
creasing pressure will cause an increasing conductance ratio. The 
conductance ratio decreases if the mobility ratio is favorable and in- 
creases if the mobility ratio is unfavorable. Ata mobility ratio of one, 
the conductance ratio will remain constant during a flood. 

Plots of conductance ratio versus areal sweep can be found in 
Appendix C. 

The ratio of the distance between injection and production wells 


ZL 
to the wellbore radius d/rw has been shown by Aronofsky and Ramey 
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to effect the conductance ratio. Although the model was constructed with 
a d/rw ratio of 320, slight differences in the well bore radius and plugging 
around the well bore will tend to change this ratio. 
The steady state flow equation for single phase flow in line drive 
24 Zoe ; 
patterns was derived by Muskat and presented by Deppe in practical 


units as 


_ 0.001538 k/uh AP 


CSRs co On eee (32) 


This equation is valid for both direct and staggered line drive well 
patterns with d/a>l. 

Using the information obtained from the model runs, the effec- 
tive d/rw ratio for the direct line drive was found to be 922. The 
effective d/rw ratio for the staggered line drive was calculated as 1110. 
These values for d/rw are somewhat smaller than those normally en- 
countered in the field. However, it is considered that reasonable 
accuracy will be obtained by using the conductance ratios obtained in this 
study without correction. 

If it is desired to correct the conductance ratios obtained in this 
study to the correct value for other d/rw ratios, using the method pre- 
sented by Aronofsky and Renee for line drives yields the following re- 


lationship: 








nn le 
' 


20 


, (rw) 
0. 682(d/a),-0. 798+In(d/ rw), +0. 434 In = 
ie = (33) 
(rw) ¢ 


where the subscripts f and m denote field and model respectively. 


Fractional Flow 


Up to the time of breakthrough, the fractional flow from the 
unswept region (oil cut) is one. Also up tothe time of breakthrough, 
a plot of areal sweep versus displaceable volumes injected yields a 
straight line at a slope of one. It can also be seen that a plot of areal 
Sweep versus elapsed time, with injection at a constant rate, will 
yield a straight line prior to breakthrough. After breakthrough, the 
areal sweep curve deviates from a straight line. 


Pega have shown that the fractional flow 


Previous publications 
from the producing well can be obtained from the slope of the areal 
sweep efficiency curve versus displaceable volumes injected. Since 
elapsed time and displaceable volumes are directly proportional, the 
fractional flow from the unswept region can be obtained by the ratio 
of the slope of the areal sweep efficiency curve at the time of interest 
to the slope of the straight line portion. Tables of fractional flow are 


given with the basic data plots in Appendix B and are plotted versus 


areal sweep efficiency in Appendix C. 





INTERPRETATION OF RESULTS 


Tables of basic data from the model runs are presented in 
Appendix B. Appendix B also includes curves of areal sweep efficiency 
versus elapsed time. Each figure also includes tables of displaceable 
volumes injected and fractional flow data as derived from these curves. 
Plots of displaceable volumes injected, fractional flow, and conductance 
ratio versus areal sweep efficiency for the mobility ratios studied are 
provided in Appendix C. 

Results of this study were compared with those obtained by the 
X-ray shadowgraph pechnianicee and were found to have a maximum 
deviation of 7% and an average difference of less than 2%. Comparison 
of areal sweep efficiencies at breakthrough and at 95% fractional flow 
from the swept region as obtained by Dyes, etal” and as obtained in 
this study are listed in Table I for the direct line drive and in Table IT 
for the staggered line drive. 

In order to make the results of this study useful for mobility 
ratios other than those specifically used here, Figures 10 through 15 
indicate the effect of mobility ratio on performance. 

In comparing the results obtained for the direct line drive and 
the staggered line drive, Figure 16 shows that areal sweep efficiencies 
are Significantly higher for the staggered line drive at breakthrough. 
However, at a fractional flow from the swept region of 95 per cent 
(assumed economic limit) the areal sweep efficiencies are similar for 
the two patterns at abandonment up to a mobility ratio of approximately 
10. At mobility ratios above 10, previous studies have shown that the 
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TAB |! 


COMPARISON OF AREAL SWEEP EFFICIENCIES 
FOR DIRECT LINE DRIVE 


Breakthrough 95% Fractional Water Flow 
Sweep Efficienc Sweep Efficiency 
Mobility 13 13 
Ratio Dyes, et al This Study | Dyes, et al This Study 
0.106 WS) 2 0.85 On 275 Ore 
0.274 Ouae 0. 76 Vee) Wess Tl 
1.0 0s>> 0.55 0.98 0.96 
313 0.46 0. 46 Opals 0.95 
Van. 1 * 0.41 * Oru 


*Not available 
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TABLE 11 


COMPARISON OF AREAL SWEEP EFFICIENCIES 
FOR STAGGERED LINE DRIVE 










Breakthrough 95% Fractional Water Flow 
Sweep Efficiency Sweep Efficiency 
Mobility 13 13 
Ratio Dyes, et al This Study 





0.106 O97 5 0799 
0.274 Oa 0.98 
1.0 0.98 Je whe 
Sig Whe 0.95 0.96 
ee 1 * 0.89 


*Not available 
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CONCLUSIONS 


The results of this study of line drive flood patterns indicate 
that mobility ratio has a definite effect on the areal sweep efficiencies 
at time of breakthrough. At abandonment conditions, the effect of 
mobility ratio on areal sweep is observable, but not to the pronounced 
extent as at breakthrough conditions. 

The conductance ratios show that the majority of the flow 
potential change takes place when the front is near the injection well 
for favorable mobility ratios and near the production well for unfavor- 
able mobility ratios. At time of breakthrough, an apparent change in 
conductance ratio is noted. 

The staggered line drive gave breakthrough sweep efficiencies 
from 4 to 14 per cent higher than those obtained by the direct line drive. 
At abandonment conditions, only small differences existed in areal 
sweep efficiencies for the direct and staggered line drive. However, 
for very poor mobility ratios, over one additional displaceable volume 
of injected fluid was required to reach abandonment conditions with the 
direct line drive. 

The method described for predicting performance of stratified 
reservoirs can be used for any well pattern for which conductance 
ratios and areal sweep efficiencies are available. This data is avail- 
able for the Recoemer. the inverted skewed Beene and the 
inverted nine-s pot’ * patterns. This study allows similar predictions 


for the direct and staggered line drive patterns. 


oy | 





The effect of a free gas saturation on areal sweep efficiency 
and conductance ratio has not been sufficiently studied. Accurate 
performance prediction for stratified reservoirs requires additional 


study on the effect of free gas saturations. 
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NOMENCLATURE 


well spacing 

constant 

row spacing of wells 

fractional area swept 

fractional flow from swept region 
fractional water flow 

height or thickness 

permeability 

length 

partial length 

mobility ratio 

fractional oil produced 

pressure drop 

total flow 

water flow 

frontal radius 

radius of radial flow region around injection well 
radius of radial flow region around production well 
radius of well bore 

saturation of connate water 
saturation of residual oil 

time 

bulk volume 


displaceable volume - Vb 6 (1 - Scw = Sro) 


a3 
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Vv velocity 

Wi water injected 

w width 

a angle open to flow into wellbore from swept region 
Y conductance ratio 

$ porosity 

nN mobility - k/u 

U viscosity 

Subscripts 

A radial flow area around injection well 
B linear flow area between radial flow areas 
C radial flow area around production well 
lesb breakthrough 

f field 

i initial 

j any stratum 

m model 

O oil 

S swept region 

t time 

u unswept region 


WwW water 
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APPENDIX A 


FLOOD FRONT TRACINGS 
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Direct Line Drive 


Vie ee Pe EE a i 


Staggered Line Drive 


FIGURE 17. FLOOD FRONT TRACINGS, M = 0.106 
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Vd 


Direct Line Drive 


Vd 


Staggered Line Drive 


FIGURE 18. FLOOD FRONT TRACINGS, M = 0.274 
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Direct Line Drive 





Staggered Line Drive 


FIGURE 19. FLOOD FRONT TRACINGS, M = 1 
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Direct Line Drive 





Staggered Line Drive 


FIGURE 20. FLOOD FRONT TRACINGS, M = 3.13 





63 


Vd 


el oo . 40 . 60 





Direct Line Drive 





Staggered Line Drive 


FIGURE 21. FLOOD FRONT TRACINGS, M = 12.1 
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TABLES AND FIGURES OF BASIC DATA 
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TABLE Il 


DIRECT LINE DRIVE 





Note: Areal sweep data is invalid. 


Run Initial Pressure Mobility Ratio 
¢ ry. 0.106 
Elapsed Corrected Fractional Pressure Conductance 
Time Time Area Ratio 
Swept 

0 = 02025 8.6 0.198 

1 ~ 0.057 9.78 0.174 

(ae O - 0.094 10.1 0.168 

STU - G. 135 WO 4 0.166 

520 = 0.181 ro. 3 0.165 
O27 - 0.672 Llc Orr 5 2 
na - 0.741 11.6 0.146 
24.27 - 0.809 4S, 0.132 
a. 27 - 0.855 136 On Z> 
A. ot ~ UPS 14.5 0.117 
Sra A - 0.937 14.9 0.114 
37.27 - 0.954 Roe, Ov) 
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TABLE IV 


DIRECT LINE DRIVE 


Run Initial Pressure Mobility Ratio 
9 2.2 0.106 
Elapsed Corrected Fractional Pressure Conductance 
Time Time Area Ratio 
Swept 
0 = 0.039 = = 
1,00 - 0.104 10.4 0.216 
Zi 16 - 0.143 10.6 0.208 
16 1h) = 0.340 10.8 0.204 
fe, 30 ~ 0.479 ae) 0.202 
20.08 “ 0.712 11.4 0.193 
23.41 ~ 0, 844 UN , 0.185 
24.41 ~ 0.870 12-6 0.175 
25.41 = - | ee 0.163 
26.41 - 0.926 14,2 O15 
27.41 - - 14.7 0.150 
29.08 - 0.949 Lon 0.144 
aD 1 > 0.974 15.4 0.143 
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TABLE V 


DIRECT LINE DRIVE 





Run Initial Pressure Mobility Ratio 
Zz bm5 0.274 
Elapsed Corrected Fractional Pressure Conductance 
Time Time Area Ratio 
Swept 

0 - 0.107 FOr? 0.607 

3333 - OF 277 lie 0.586 

6.67 = 0.453 1 brs 0. 505 
10.00 - 0.697 Vane 0.520 
i325 - 0. 860 14.1 0.461 
16.58 ~ 02925 14.9 0.436 
23.08 7 0.945 L520 0.433 
Zions 3 - - ey, 0.425 
29.67 - 0,970 15.4 0.422 
Note: Rerun as run 22. 
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TABLE VI 


DIRECT LINE DRIVE 


Run Initial Pressure Mobility Ratio 
9 4.8 0.274 
Elapsed Corrected Fractional Pressure Conductance 
Time Time Area Ratio 





0 0 0.194 emo 0.658 
1,0 a 0.288 to 0.658 
2.0 30 0.372 133 0.658 
3.0 Dat 0.467 te3 0.658 
4.0 t0 Ono o8 y ore. 0.658 
5.0 aa 0.665 Tes 0.658 
fa)0 13. 0.848 og 0.565 
8.0 Nya 7 F0 02533 
9.0 fecal Orato) Te Orne OF exe 74 
10.0 ees - yee. Oe203 


11.0 AU, 0.930 
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S86°0 S6L°I 026°0 0°92 
€96°0 2L299°T 296°0 0°r2 
€€6°0 O2S°T 196°0 0°22 
L68°0 28€°I 1596°0 0°02 
1€8°0 €F2'T 9€6°0 0°ST 
O€2°0 SOT‘*I O0T6°0 0°9T 
845°0 296°0 028°0 O°FI 
SZ2€°0 628°0 O18°0 Oct 
000°0 092°0 092°0 O°1l 
ay PA Sa SUIT] 
xO 
Le) ae . 


cl 


GWIL GASdVTd sustoa LqdaMS VAUV -- AAIYGC ANIT LOAYIAG 


a nals | 


cUnNY Oo 


"€2 AUN 





TABLE VII 


DIRECT LINE DRIVE 








Run Mobility Ratio 
] 1.0 
Elapsed Corrected Fractional Pressure Conductance 
Time ‘Gare Area Ratio 
Swept 
0 - 0.0944 
lee - 0, 240 C : 
2,0 ~ 0.385 O 
a0 - OFS19 N 
4.0 - 0.639 S ° 
550 ~ 0. 726 T 
6.33 7 On77t A 
i 33 ~ 0. 840 N 
10.58 ~ 0.939 E 
LoS - 0.973 Lee 


Note: Rerun as run Zl, 
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TABLE VIII 


DIRECT LINE DRIVE 


Run Mobility Ratio 
6 1,0 
Elapsed Corrected Fractional Pressure Conductance 
Time Time Area Ratio 
Swe pt 

0 0 Oued? 1 1.0 

1.0 0,84 0. 356 C 

220 1.68 0. 450 O 

3.0 2.52 0.546 N le 

4.0 3.36 0.605 S le 

5». 0 4.21 ON00 / a 1 

6.0 5,05 0.724 A i 

8.0 6.72 0.790 N 1. 
10.0 8.42 0.845 ae i 
me. O97 | 0, 0.934 i 
elqisky, 15.60 0.960 ie 
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we “ae 
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TABLE IX: 


DIRECT LINE DRIVE 


Run Initial Pressure Mobility Ratio 
3 s.3 3u13 
Elapsed Corrected Fractional Pressure Conductance 
Time Time Area Ratio 
Swept 
0 ~ 0.099 = = 
170 - 0,210 - - 
2.0 ~ Ove? eae! 1, 46 
3.0 ~ 0.432 Oe 1, 46 
4.0 - 0.482 Si. 1,64 
5.10 ~ 0,530 ice, 1e3s 
lore 78 - 0.654 te0 Lace 
17.40 - 0.750 1.0 i 70 
Zl, 25 - 0.798 TO Lege iG, 
28.00 - 0.862 7.0 ie 0 
34.00 ~ 0. 896 tao 1.90 
35.00 - Omens 420 1.90 


Note: Areal sweep data is invalid. 
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TABLE X 


DIRECT LINE DRIVE 


Run Initial Pressure Mobility Ratio 
8 12.4 oa eS 
Elapsed Corrected Fractional Pressure Conductance 
Time Time Area Ratio 
Swept 

0 - 0.051 -~ - 
10 - Oreo 9.4 leo Z 
ZO - Ueo50 G25 1.46 
2.0 - = Sree Hee oz, 
4.0 : 0.580 7.0 lea? 
5.0 - 0.642 655 et 
6.0 ~ - 6.3 ee el 
0 - 0.720 6.2 2.00 
12.08 - 0.829 6.0 ZU 
lezo8 - 0.920 6.0 2.07 
22.42 - 0.952 6.0 2.07 
25.50 - 0.965 6.0 2.0 / 
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TABLE XI 


DIRECT LINE DRIVE 


Run Initial Pressure Mobility Ratio 
5 14.7 2 
Elapsed Corrected Fractional Pressure Conductance 
Time Time Area Ratio 
Swept 

0 ~ 0.065 = - 

1 20 - 0.101] ~ 7 

2.0 - Omr5 7 On 7 ael9 

3.0 - 0.207 oO. 1 rap, 

4.0 - 0. 259 Oe? Fe li 

6.0 - 0.354 ery Fay NN 

a0 - 0.410 6.6 2.23 

8.0 ~ - oe, 2.49 

o,.0 - 0.496 B10 2.94 
10.0 7 ~ 4.4 3.34 
11.0 - 0525 a5 4.20 
Vas0 - ~ Zao Sol? 
13:0 - 0.545 Zac 5), Ae. 
SRro 2 ~ 0.586 Za 5.45 
22093 = - 2.6 505 
Ag eS) ~ 0.616 Ze 6.40 
39.67 - 0.670 Zi3 6.40 
5on 03 - 0. 709 Zo 6. 40 
74.25 “ 0.742 Ze 6.40 
oi - ORT Tre: Ze 6.40 
Nay ao ie - 0.812 25 6.40 
Note: Areal sweep data is invalid. 
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fieJeliid, ><ige 


DIRECT LINE DRIVE 


Run Initial Pressure Mobility Ratio 
10 1.0 yee 
Elapsed Corrected Fractional Pressure Conductance 
Time Time Area Ratio 
Swept 

0 0.038 ~ 
10 0.107 one 1,65 
Ze 0 Oeets : ae 5 
3.0 0.247 1. 70 
4.0 02323 : Leaf 
55.0 0.401 li) 
6.0 = : 2.14 
7.0 0.507 4. 3. om 
8.0 ~ 376 4.17 
9.0 O2551 ‘ 4.41 
10.0 - 4 4.4] 
12,0 0.654 : 5.00 
24.59 0.706 50 
28.34 ~ : 5,750 
34.92 0. 740 : Scie 
41.50 0. 169 sie Sil, 
ala 0.807 : 5.36 
62. 25 0.833 2.8 5280 
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TABLE XIII 


STAGGERED LINE DRIVE 


Run Initial Pressure Mobility Ratio 
14 Zen 0.106 
Elapsed Corrected Fractional Pressure Conductance 
Time Time Area Ratio 
Swept 
0 - ae dk sf 0.230 
6.08 - 0.282 im 9S 0.229 
13.08 = 0.482 HZ. 0 0.225 
16.66 ~ Om59S ay i 0.223 
23.16 - 0. 768 2, 5 0.26 
Zo. 16 - 0.918 lh arg 0.197 
oy ..16 - ~ 14.2 0.190 
28.16 ~ O95 5 14.7 0.184 
29.16 - 02967 ie 0 0.180 
30.16 = - Ey A 4 0.178 
Sal 6 = - 15.4 O75 
32416 ~ 0.978 1.6 Onl 2 
BY 16 - 0.989 Woe F Ot 
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TABLE XIV 


STAGGERED LINE DRIVE 


Run Initial Pressure Mobility Ratio 
le, (a) 0.106 
Elapsed Corrected Fractional Pressure Conductance 
ame Time Area Ratio 
Swept 

0 0 0.0086 9.6 0.284 

e8 0.505 0.032 Ze 0.231 

DcO (hg en 0.074 12.2 0.224 
1 fac Mey Srl 0.174 LZ. 0.224 
aero 2 7.94 0.279 yaa, 0.224 
Zones 10.80 0.384 b2E5 0.224 
32.05 13265 0.503 Vee 0 0.224 
38.62 16.45 0.630 EZa5 0.224 
45.20 19. 20 re ar, PZ 5 0.224 
51.78 22.00 0.822 BARES, 0.224 
So) lA Zo, 50 0.867 13720 O52l5 
58.45 24.90 0.900 14.0 0.200 
61.70 26.30 0,932 10 0.187 
64.95 27.60 0.948 [5.6 Ont 
74. 70 31.80 On > 1/5) as: OLE Z 
84.37 355.90 Ong 90 16.8 0.170 
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TABLE XV 


STAGGERED LINE DRIVE 





Run Initial Pressure Mobility Ratio 
lh oes, 0.274 
Elapsed Corrected Fractional Pressure Conductance 
Time Time Area Ratio 
Swept 
0 - 0.066 - ~ 
a. 0 - OLS ~ - 
4.0 ~ OnZ75 - - 
a0 - 0.347 ~ ~ 
6.0 0.396 Oy 0, Sas 
Gao - 0. 461 10.7 0.514 
0 = 0.584 WOR t 0.514 
11,0 - 0. 76 LO 0.514 
or 0 - 0.842 1 ae 0.486 
14.0 - OF 725 Lao 0.466 
po R0 - 0.950 12.0 0.458 
ii. oO a 0.971 12.2 0.451 
ZO - 0.980 NAR 0.45] 
250 - 0.988 12.2 0,451 
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TABLE XVI 


STAGGERED LINE DRIVE 


Run Initial Pressure Mobility Ratio 
17 oral 0.274 
Elapsed Corrected Fractional Pressure Conductance 
Time Time Area Ratio 
Swept 

0 0 0.414 oll 0.560 
ia. 0 0.934 OF 371 Sho & 0.537 
2.0 lene: 0.143 oie 2) 0), Sy Sig 
7.0 280 0.196 9.6 Oma 51 
4.0 3.74 Os Aas 9.6 0. el 
6.0 5.60 Oreoy 2 9.6 O@>5 | 
8.0 7.47 0.481 9.6 One 
7.0 9.34 Ve eheare 9.6 0.531 
re. 0 Dice 0. 703 Shee fl 0.526 
MSA ba.20 - 26 Oma20 
14.0 P3520 Oe EES, POR 02505 
15.0 14.00 0.884 10.6 0.481 
V6. 0 14.94 - Li 0 0.464 
i. 0 15.88 Oe Bl3 0.451 
EO 15 0.976 Na res 0.443 
2020 18.68 - ie Bas: 0.440 
ZY 20 19.60 0.987 iy 0.436 
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TABLE XVII 


STAGGERED LINE DRIVE 


Run Mobility Ratio 
11 Lesa) 
Time Time Area Ratio 
Swept 

0 0 0.083 C Lae 
0) Peo Oa Tay O 0 
Z. 0 Ze 0.473 N Fs 0 
S10 het a) 0.685 S a 
4.0 4.20 0-625 a L508 
5) 510) As, 0.904 A 1.0 
6.0 6.30 Oe N ba8 
3) 0 8.45 0.987 T ra) 
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TABLE XVIII 


STAGGERED LINE DRIVE 


Run Mobility Ratio 
16 1.0 
Elapsed Corrected Fractional Pressure Conductance 
Time Time Area Ratio 
Swept 
= O; 113 C Le 
1.0 ~ Oez51 O i 
Z.0 ~ 0.450 N ie 
3.0 = 0.611 S i 
4.0 - a af ae) a ie 
5.0 - 0.876 A ys 
6.0 - 0.940 N I 
7.0 - 0.967 S i 
8.0 - 0.992 T LNs 
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TABLE XIX 


STAGGERED LINE DRIVE 


Run Initial Pressure Mobility Ratio 
13 14.7 a 
Elapsed Corrected Fractional Pressure Conductance 
Time Time Area Ratio 
Swept 
0 - 0.098 Se 62 
es - OA, 8.9 1265 
2:0 - 0. 369 8.9 16> 
aU - 0.509 8.8 leo 
4.0 ~ 0.636 8.8 1.67 
5.0 - On139 Sel 1.81 
6.0 ~ 0.778 liars Z.01 
7.0 ~ - gira fae 0) Ty 
8.0 - 0.832 Geo Zao 
10.0 - 0.873 C0 Zu 
13.0 - ee 7.0 Fase MG, 
5), O - 0.936 CO Caeu0 
17.0 ~ 0.969 7.0 2.10 
19.20 - 0.988 rma 210 





Run 
18 


Elapsed 


Time 


0 
0 
5.0 
4.0 
on © 
6.0 
aw 
8.0 
920 
10.0 
om. 0 
12.0 
rs. 0 
15.0 
Pond 
20.0 
Z 32.0 


TABLE XX 


STAGGERED LINE DRIVE 


Corrected 
Time 


0.728 


1, 46 
Qe) 
2.92 
3.64 
Ai, 
5.10 
Sore, 
6.56 
UG TAS: 
S202 
Sint 
9. 48 
10.92 
13.0 
14, 60 
16. 80 


Initial Pressure 
14.2 


Fractional 


Area 
Swept 


0.098 
05191 
0.287 
OS Wat 
0.465 
O55 5 
BN i6i74 
0.694 
O51 59 
0,802 


02657 
0. 866 
O3595 
0.920 
0.956 
ee ee 


Mobility Ratio 
Se 


Pressure 


LOS. 
TO 
10.4 
Gee 
9a9 
ee, 
Ore 
ae 
ee 
9.0 
ae 
7.90 
7 aw 


Ratio 


15 
1535 
Lest 
1.41 
1.43 
1.49 
Wa Sys. 
1/56 
Leos 
[253 
Mr Sie: 
iS 
1,56 


90 


Conductance 





eh 





cl°c=W 
ANIL CGASdVTA snst94 LGAMS VAUV -- AAINC ANIT CAUYSBDOVLS ‘o¢ TUNOM 
SULT, 
02 91 ZI g v 


stTuny x 
ee | Stsl 1G) 





OZ 
TABLE XXI 


STAGGERED LINE DRIVE 


Run Initial Pressure Mobility Ratio 
Is, 16.0 Ae | 
Elapsed Corrected Fractional Pressure Conductance 
Time Time Area Ratio 
Swept 

0 - 0.014 - - 

1.0 - 0.072 8.6 1,86 

2.0 - 0.134 8.2 he, 

3.0 = - A, 2.05 

4.0 - 0.254 fe Zoe 

22,0 - 0.315 7.4 ZaN6 

7.0 - 0.424 Maas a 29 

8.0 - 0.477 7.0 Zae9 

20 - . 62/3 2.54 
10.0 - 0.584 4.8 3.33 
12.0 - 0.644 4.8 3233 
13.0 - - 4.2 ere) 
14.0 ~ ~ 33 4.21 
16.0 - 0.715 3.4 4.71 
18.0 - 0.749 3.4 4.71 
24.91 ~ 0.785 3.4 4.71 
42.65 - 0.877 3.4 4.71 
55.08 - 0.915 3.4 4.71 
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TABLE XXII 


STAGGERED LINE DRIVE 





Run Initial Pressure Mobility Ratio 
Z0 LO. > I ge | 
Elapsed Corrected Fractional Pressure Conductance 
ite Time Area Ratio 
Swept 

0 0 0.038 ie UP is, 

Fe O OQ 525 Op 04 oe 2.01 

4.33 rapa | 0.204 ; ele 

m. 58 Be 0.294 a8 Zn 
Or. 9 1 6.79 0.385 (exe: Paes Ml 
14.16 7.42 0.468 7.8 Z.11 
i. 5O J 0.553 7.8 Cee 
18.50 9.68 0.57) Ure age lal 
ZUR O 10.7 - led 2.20 
Zr. 50 2 0.648 TA Vag SS, 
Za. 30 Ee! ~ (20 2530 
25.00 13.1 0.693 6.8 2.43 
Boe co 14.8 - Gee. 2.66 
31.66 16.6 0.747 shes) 3 10 
SN AS) 20.5 Oar! Dian oa) 
48.08 Als Fe | 0.807 5.5 320 
3) T) Age Re, SOR 0.827 as, 3.0 
67.41 35,.2 0.853 ores. S70) 
76.08 39.8 0.872 Bien 3.0 
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APPENDIX C 


FIGURES OF INTERPRETED DATA 
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FIGURE 32. DIRECT LINE DRIVE -- MOBILITY RATIOS 
for DISPLACEABLE VOLUMES INJECTED versus 
AREA SWEPT 
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FIGURE 33. 
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DIRECT LINE DRIVE -- MOBILITY RATIOS 
for FRACTIONAL FLOW versus AREA SWEPT 
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FIGURE 34. DIRECT LINE DRIVE -- MOBILITY RATIOS 
for CONDUCTANCE RATIO versus AREA SWEPT 
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FIGURE 35. STAGGERED LINE DRIVE -- MOBILITY RATIOS 
for DISPLACEABLE VOLUME INJECTED versus 
AREA SWEPT 
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FIGURE 36. STAGGERED LINE DRIVE -- MOBILITY RATIOS 
for FRACTIONAL FLOW versus AREA SWEPT 
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FIGURE 37. STAGGERED LINE DRIVE -- MOBILITY RATIOS 
for CONDUCTANCE RATIO versus AREA SWEPT 





APPENDIX D 


COMPUTED CONDUCTANCE RATIOS 


is. 


nu se =e aoe +e 
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Os 


PROGRAM CONRA 
DIMENSION RM(10) »DFS(20) »ES(10920) 

READ ls JEs(RM(J)»s J=19JE) 

READ ls KEs (DFS(K)»s K=1sKE) 

FORMAT (13911F7.0) 

Per 10 J=1.JE 

FORMAT (11F7.0) 

READ 29 (ES(JsK)»s K=1 KE) 

RM IS THE MOBILITY RATIO AND JE IS THE NUMBER OF MOBILITY RATIOS 

WHICH ARE TO BE STUDIED. DFS IS FRACTIONAL FLOW FROM THE SWEPT REGION 
AND KE IS THE NUMBER OF FRACTIONAL FLOW VALUES TO BE READ IN. 

ES VALUES ARE AREAL SWEEP EFFICIENCIES CORRESPONDING TO THE DFS VALUES. 
SET INITIAL VALUES 


D 


1.0 

A =D 

REI = 0.25*D 

RW = D/320.0 

RLW= O0«5*(D/A ~—- 14617 + 06637*LOGF(A/REI) ) 
ESA = 3-1416*REIT**2/(2.0*D*A) 

Domo" J=15JE 


PRINT 3sRM(J) 


Z3O0FORMAT (1H1940Xs16HMOBILITY RATIO =F5e2//10Xs19H CONDUCTIVITY RATI 


10910Xs14HFRACTION SWEPT 913X92HFS//23Xs7H1 00000 219X 91HO 219X 91 HO) 


Eobd = ES( Js) 


REP = SQRTF(2.0*D*A*(1.0-ESBT)/3.1416) 
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DES = 0-01 
ES(Jo1) = DES 
DO 70 K=15KE 
80 IF (ES(JsK) - ESA) 20520930 
130 IF (ES(JsK) - ESBT) 40550,50 
20 DLW = 020 
RF = SQRTF(2.0*D*A¥ES (JK) /3.1416) 
GO TO 60 


40 RF = REI 


DLW =RLW*(ES(J9K) - ESA)/(ESBT - ESA) 
DES = 0-05 
GO TO 60 
BO RF = REI 
DLW = RLW 
600CR = (LOGF(REI/RW) + 6628¥*RLW+ LOGF(REP/RW))/(1-¢/RM(J) *LOGF (RF /RW) 
1+ LOGF(REI/RF) + 6¢28%(1-e/RM(J) - 120) *DLW + 6628*RLW+ LOGF (REP/RW 


2)*(1.0 — DFS(K))/(1.0-DFS(K)/(RM(J)—-DFS(K)*(RM(J)-160)))) 
PRINT 45 CRsES( JK) sDFS(K) 
4 FORMAT (10X»3F20.5) 
ES(JsK) = ES(JsK) + DES 
IF (ES(JsK) - ESBT) 80570270 
70 CONTINUE 


END 


oe atest <= 








Ke 
Flow Diagram Program CONRA 






Dimension 
M, fs, Es 
Read 
M, fs, Es 


Set initial values 


d=1.0 







rw = d/ 320.0 
ay Bus 0}: o(a/a - 1.17 + 0. 637 loga/ri ) 
Es, iii ri 9 Z2da 





Do Loop 
for all 
M values 


= Es(j, 2) 


AEs = 0.01 


Es(j,1) = SEs 


Do Loop 
for all 
fs values 
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eG =e 
fs = 0 
rf =(2da es) i oe 













ra 
0 
=10.,05 


Liw _L(Es - Esa) 


w(Es, 7" Es ,) 









Y = (Inri/rw + 6.28 L/w + Inrp/rw) / 
((1/M) In rf/rw + lnri/rf + 6.28(1/M - 1) 4/w 
+ 6.28 L/w + lnrp/rw [1 - fs / (M- fs(M - 1))] 


Print 
VY, Es, fs 


Es = Es + AEs 


Oo 
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FORTRAN SYMBOL CONVERSION CHART 


FORTRAN 


A 
CR 
D 
DES 
ES 


ESA 


ESBT 
DFS 
DLW 


REI 


REP 


RF 
RLW 
RM 


RW 


Nomenclature 


a 


L/w 
ri 


rp 


rf 


L/w 


rw 


Description 
well spacing 
conductance ratio 
row spacing of wells 
increment of area swept 
fractional area swept 


radial flow region around 
injection well 


area swept at breakthrough 
fractional flow from swept region 
partial length to width ratio 


radius of radial flow region 
around injection well 


radius of radial flow region 
around production well 


frontal radius 
length to width ratio 
mobility ratio 


radius of wellbore 
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— —— Analytical 
—— Model 


FIGURE 38. DIRECT LINE DRIVE -- COMPARISON OF 
ANALYTICAL and MODEL CONDUCTANCE RATIOS 
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FIGURE 39. STAGGERED LINE DRIVE -- COMPARISON OF 
ANALYTICAL and MODEL CONDUCTANCE RATIOS 
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FIGURE 40. DIRECT LINE DRIVE -- EFFECT OF L/w RATIO 
ON CONDUCTANCE RATIO 
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EXAMPLE PREDICTION FOR STRATIFIED RESERVOIR 





APPENDIX E 


EXAMPLE PREDICTION FOR STRATIFIED RESERVOIR 


To illustrate the procedure described for stratified reservoirs, 
assume a staggered line drive unit of 40 acres well spacing composed 
of 2 strata. Under reservoir conditions, the oil viscosity is 4 cp and 
the water viscosity is 0.4cp. Assume that no gas saturation exists. 


The reservoir properties are: 


Stratum l Stratum 2 
Height 20 feet 20 feet 
Porosity 0. 30 O25 
Connate water (Scw) 0. 30 OF 25 
Residual oil (Sro) 0. 30 0.30 
KO (at Scw) 200 md 75 md 
ky (at Sro) 20 md 22 md 


The bulk volume in each stratum will be the same or (20) (80) = 


1600 acre ft. 


Vd = Vb ¢$(1-Scw-Sro) 

Vd, = (1600) (0. 30) (1-0. 30-0. 30) = 192 acre-ft 

Vd, = (1600) (0. 25) (1-0.25-0.30) = 180 acre-ft 
The mobility ratios for each stratum are 

M = (kK. /y)/(k./u,) 

M, = (20/0. 4)/(200/4) = 1 

M, = (22/0. 4)/(75/4) = 3.13 


2 





ei (2oenroe 


= 2.84 
2Vd> 


To obtain ct for each stratum from equation (25) the inverse 
ofthe conductance ratio is plotted versus displaceable volumes injected 
(from Figures 13 and 15). The resulting curves (Figure 41) are 
graphically integrated to obtain ct, and ct, . (Note that the data for 
each Figure is in the preceding Table. ) 

Cumulative water injected and the corresponding oil produced 
are obtained from Figure 13 in dimensionless terms of displaceable 
volumes injected and fractional area swept. These are then converted 
to Wi and Np and plotted for each stratum versus ct in Figure 42, 
Total water injected and oil produced can be obtained from Figure 42 
by summing the Wi and Np curves for each stratum and is plotted 
in Figure 43, This curve can be converted to an actual time relation- 
Ship by knowing the maximum injection rate and well allowables. 
Fractional flow from the swept region is obtained as the slope of the 
Np versus Wi curve and the total conductance ratio is obtained from 
equation (31). Fractional flow is plotted on Figure 43 and conductance 
ratio versus water injected for the system is plotted on Figure 44, 

The above example does not consider the presence ofa free gas 
saturation at the start of the water flood. However, it will be within 
the realm of engineering accuracy to assume that oil bank fill up occurs 
simultaneously in all strata. This will be true except in the case of 
high gas saturations because of the unusually favorable mobility ratio 


during oil bank fill up. The volume occupied by free gas prior to oil 
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bank fill up can then be added to the cumulative water injected and 


volume swept to obtain satisfactory performance predictions when there 


is free gas present at the beginning of the flood. 





Wi/Vd 


fraction 


0.1 
RS) 
0.065 
1.0 
| ee) 
2.0 
faye oNs) 


TABLE XXIII 


acre-ft 


18.0 
705.0 
6 es: 
180.0 
270.0 
360.0 
450.0 


Y2 


1.0 

1.48 
ae ¢. 
152 
Meads, 
Nees, 
1.84 
WE o)8 


l/y, 


1G 

Ome 
0075 
Ovosy 
0.501 
0.546 
0,544 
0.540 
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TABLE XXIV 


Stratum one Stratum two 

Wi ct) 2 
Wi (l/y) dWi 2.84 (1/y) dWi 

0.0 0.0 0.0 0.0 0.0 

100.0 100,0 50.0 35. 38 100.5 

200.0 200.0 100.0 68.83 IW) 5: 

300.0 300.0 150.0 100. 88 285.0 

400.0 400.0 200.0 132. 88 377.0 

300.0 188.18 55 52.0 

400.0 242.48 688.0 

500.0 296. 38 843.0 

Wi/Vd Np/Vd Wi Np Wi/Vd Np/Vd Wi Np 

fraction fraction acre-ft acre-ft |} fraction fraction acre-ft acre-ft 
0.769 0.769 147.5 147.5 0.665 0.665 119:5 119.5 
0.80 C775 [S35 P522-5 0.80 0.743 143.8 133.5 
e270 0.90 192.0 173.0 1.0 0.803 180.0 144.2 
lee Z OF 752 Z50-0 182.8 es, 0.883 269.0 L56may 
1,4 02935 ZoonO 189.0 ZO 0.930 359,0 L670 
Cao 0.961 449.0 173.0 
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ct 


100 
145 
200 
225 
250 
300 
350 
400 
500 
600 
700 


Wi 


100 
145 
200 
225 
250 
300 
B20 
400 
500 
600 
700 


TABLE XXV 


Np 


Np 





Z ] 2 
50 100 50 
75 145 15 
105 ay ge 105 
120 180 120 
132 185 2a 
160 190 138 
rey 190 146 
216 190 L152 
276 190 160 
344 180 166 
410 190 171 
Wi fo fw 
Z20 1.000 0.000 
300 05570 0.410 
340 0.550 0. 450 
350 0.400 0.600 
400 0.245 Oe 15> 
500 Oeee? 0.873 
600 0,054 0.946 
800 0.038 0.962 
1000 07031 0.969 





Wi 


150 
220 
305 
345 
382 
460 
537 
616 
CaS 
944 
0 


Np 


150 
220 
278 
300 
Sah2 
328 
336 
342 
350 
356 
361 


by 


) 
4 

} 
e 
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ct 


50 
100 
150 
200 
300 
400 
500 
600 
700 


11 teva .e) 
1+ peat 
pee] 
Wi, l/y> 
25 04075 
50 0.672 
78 0.670 
105 0.660 
160 0.600 
216 On 552 
Zio 0.550 
344 0.542 
410 0.540 


TABLE XAXVI 


k,Vd 


Oe Sen. 
geo (es 7A ad 
Y> Wi 
1.480 75 
1.490 150 
1.492 225 
1.518 305 
1,669 460 
1.812 616 
1.820 776 
1.845 944 
1.850 1110 


. 126 
.128 
.128 
134 
« ta 
<c 80 
214 
. 220 
me 4 ih 


iLvaglk 
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